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The rapid transition between spherical and 7-soft shapes in Ba and Xe nuclei in the mass region 
A > 130 is analyzed using excitation spectra and collective wave functions obtained by diagonal- 
ization of a five-dimensional Hamiltonian for quadrupole vibrational and rotational degrees of free- 
dom, with parameters determined by constrained self-consistent relativistic mean-field calculations 
for triaxial shapes. The results reproduce the characteristic evolution of excitation spectra and E2 
transition probabilities and, in general, a good agreement with available data is obtained. The cal- 
culated spectra display fingerprints of a second-order shape phase transition that can approximately 
be described by analytic solutions corresponding to the E(5) dynamical symmetry. 
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I. INTRODUCTION 

Atomic nuclei display a variety of different equilibrium shapes - spherical, axially deformed, shapes that are soft 
with respect to triaxial deformations. In most cases the transitions between different shapes in isotopic or isotonic 
sequences are gradual, and reflect the underlying modifications of single-nucleon shell structure and interactions 
between valence nucleons. In a number of cases, however, with the addition or subtraction of only few nucleons one 
finds experimental signatures of abrupt changes in observables that characterize ground-state nuclear shapes. In the 
last decade the concept of quantum phase transitions (QPT) has successfully been applied and investigated, both 
experimentally and theoretically, in equilibrium shape changes of nuclei [1-6]. 

The two best studied classes of nuclear shape phase transitions correspond to a second-order QPT between spherical 
and 7-soft shapes [7, 8], and a first-order QPT between spherical and axially deformed shapes [9, 10]. Most theoretical 
studies of QPT in nuclei have been based on model specific Hamiltonians that by construction describe shape changes, 
i.e. a phase transition is accessed by variation of a control parameter explicitly built into the Hamiltonian. For instance, 
a QPT can be described in the geometric framework in terms of a Bohr Hamiltonian for shape variables, and related to 
the concept of critical symmetries that provide parameter independent predictions for excitation spectra and electric 
quadrupole (E2) transition rates for nuclei at the phase transition point. Analytic solutions of the eigenvalue problem 
at the critical point are associated with zeros of special functions. Alternatively, in the algebraic approach different 
shapes coincide with particular dynamic symmetries of some algebraic structure, and a QPT may occur when these 
symmetries are broken in a specific way. An example of the latter approach is the framework of the intcracting-boson 
model (IBM) [11] and its various extensions. 

More recently a number of studies appeared that attempt a fully microscopic description of shape QPT starting from 
nucleonic degrees of freedom [12-20]. An approach in terms of explicit nuclconic degrees of freedom is particularly 
important because the physical control parameter in a nuclear QPT is the actual number of nucleons. In a series 
of studies [15, 19, 20] we investigated nuclear QPT in the region Z = 60, 62, 64 with N sa 90, using a microscopic 
approach based on constrained self-consistent relativistic mean-field calculations of potential energy surfaces. Since 
a quantitative analysis of QPT must go beyond a simple mean-field level, i.e. one must be able to calculate ratios 
of excitation energies and electromagnetic transition rates, in Ref. [15] the generator coordinate method (GCM) 
was used to perform configuration mixing of angular-momentum and particle-number projected relativistic wave 
functions restricted to axial symmetry. This approach was extended in [19, 20], where collective excitation spectra 
and transition probabilities were calculated starting from a five-dimensional Hamiltonian for quadrupole vibrational 
and rotational degrees of freedom, with parameters determined by constrained mean-field calculations for triaxial 
shapes, i.e. including both /3 and 7 deformations. The results reproduced available data, and have shown that there 
is an abrupt change of structure at = 90 that can approximately be characterized by the X(5) analytic solution [9] 
at the critical point of the first-order quantum phase transition between spherical and axially deformed shapes. As a 
function of the physical control parameter - the number of nucleons, energy gaps between the ground state and the 
excited vibrational states with zero angular momentum, isomer shifts, and the monopole transition strength, exhibit 
sharp discontinuities at neutron number N ~ 90, characteristic of a first-order quantum phase transition. 

Are the remarkable results for X(5)-like shape phase transitions in the mass region A w 150 [15, 19, 20] somewhat 
accidental? Can the same universal energy density functional describe other types of shape phase transitions, and in 
different regions of the chart of nuclides? To answer these questions we apply the same model to the study of 7-soft 
shapes in Ba and Xe nuclei in the mass A > 130 region, where evidence has been reported for a second-order QPT 
between spherical and 7-soft shapes. This is a phase transition in one degree of freedom - the axial deformation /3 
and, in the IBM language, represents a transition between the U(5) and 0(6) dynamical symmetries in the limit of 
large boson number. The critical point of phase transition can also be related to a dynamical symmetry: in this case 
E(5) [7], and the experimental realization of this critical-point symmetry was first identified in the spectrum of ^"^^Ba 
[8]. 

The model is based on constrained self-consistent relativistic mean-field calculations for triaxial shapes, i.e. in- 
cluding both (3 and 7 deformations. The resulting self-consistent solutions: single-particle wave functions, occupation 
probabilities, and quasiparticle energies that correspond to each point on the binding energy surface, arc used to 
calculate the parameters that determine the collective Hamiltonian: three mass parameters, three moments of inertia, 
and the zero-point energy corrections, as functions of the deformations /3 and 7 [21]. The diagonalization of the 
Hamiltonian yields the excitation energies and collective wave functions that are used to calculate observables. No 
attempt is made to tune the model or the parameters to this particular type of shape transition or mass region. 
As in our previous studies of X(5)-like transitions [15, 19, 20], the relativistic functional PC-Fl (point-coupling La- 
grangian) [22] is used in the particle-hole channel, and a density-independent (5-force is the effective interaction in the 
particle-particle channel, with pairing correlations treated in the BCS approximation. 
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II. SHAPE TRANSITION IN Ba AND Xe ISOTOPES 



Our analysis starts with the self-consistent relativistic mean-field plus BCS calculation of quadrupole energy surfaces 
of Ba and Xe isotopes. The map of the energy surface as function of the quadrupole deformation is obtained by 
imposing constraints on the axial and triaxial mass quadrupole moments. The method of quadratic constraints uses 
an unrestricted variation of the function 



(H) + J2 ((02^) - 92 



(1) 



where (H) is the total energy, and {Q2fj.) denotes the expectation value of the mass quadrupole operators: 
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q2^ is the constrained value of the multipolc moment, and C^p, the corresponding stiffness constant [23]. 

In Fig. 1 we display the RMF+BCS triaxial quadrupole binding energy maps of the even-even ^'^^~^'^°Ba in the /3 — 7 
plane (0 < 7 < 60°). All energies are normalized with respect to the binding energy of the absolute minimum, and 
the contours join points on the surface with the same energy (in MeV). For each nucleus in the corresponding inset 
we plot the axial projection of the binding energy for oblate (negative (3) and prolate (positive /?) deformations. The 
corresponding energy maps of the even-even ^■^^~^^^Xe isotopes are shown in Fig. 2. Both for Ba and Xe, these plots 
illustrate the rapid transition from spherical shapes near the N = 82 closed-shell, to 7-soft energy surfaces for lighter 
isotopes. Starting from almost perfect spherical shapes in ^'^^Ba and ^"^^Xe, the decrease in neutron number induces 
quadrupole deformations. In the axial plots we notice that the potentials become less stiff in /3 and that two minima 
develop, one prolate and one oblate, at almost the same energy, and separated by a low spherical barrier. However, 
the 3-dimensional energy maps show that the spherical barriers are concentrated around /3 « and, therefore, rather 
than two separate minima, the potentials display continuous 7-soft minima that extend from prolate to oblate shapes. 

Of particular interest in the present analysis are the nuclei which have been identified as possible candidates for a 
shape phase transition that can be characterized by the E(5) dynamical symmetry [7]. The experimental realization 
of this critical-point symmetry, associated with a second order quantum phase transition between spherical and 7-soft 
potential shapes, was first identified in ^'^^Ba [8]. E(5) is the symmetry of a 5-dimcnsional (intrinsic variables /3 and 7 
and the three Euler angles) infinite well in the axial deformation variable (3 ( V{/3) — for \(3\ < /3w, and V{P) = 00 
for |/3| > Pw), and the potential is completely 7-independent. The microscopic binding energy curve i?(/3) of ^^^^Ba 
(Fig. 1) displays a shape that is almost symmetric with respect to /3 = 0. One notices a relatively flat bottom between 
/? « —0.1 and /? « 0.1 (the oblate configuration is only « 0.5 MeV above the prolate minimum), and the potential 
is rather stiff for \(3\ > 0.15. The dependence on the triaxial deformation parameter 7 is seen in the corresponding 
3-dimensional energy map and, even more clearly, in Fig. 3 where we plot the binding energy curves as functions of 7 
for several values of the axial deformation: (3 = 0.05, 0.1, 0.15, and 0.2. In the region of the flat bottom |/3| < 0.1 the 
binding energy of -'^'^^Ba is indeed almost independent of 7, and even for somewhat larger deformations 0.1 < |/3| < 0.2, 
only a weak dependence on 7 is predicted by the calculation based on the PC-Fl functional. A very similar energy 
surface is calculated for the isotone ^■^^Xe (Figs. 2 and 4). 

In the next step the constrained self-consistent solutions of the relativistic mean-field plus BCS equations, i.e. 
the single-particle wave functions, occupation probabilities, and quasiparticle energies that correspond to each point 
on the binding energy surfaces of Figs. 1 and 2, are used to calculate the parameters that determine the collective 
Hamiltonian [21]: 



H — Tvib + Trot + V^coll 



(3) 



with the vibrational kinetic energy: 
Tvih 



and rotational kinetic energy: 



1 

1 



/? sin 37 



d I r 
d I r 



dp V w^^^'^'^dp dp\l u;^'^^^97 



3 ~ df: 
d 1 d 



sin3 J Bf}^ , 

op P \ w 



(4) 



2 ^ Ik 

k=l " 



(5) 



4 



Vcoi\ is the collective potential. Jk denotes the components of the angular momentum in the body-fixed frame of a 
nucleus, and the mass parameters S^/j, Bp^, B^^, as well as the moments of inertia 1^, depend on the quadrupole 
deformation variables /3 and 7: 



Two additional quantities that appear in the expression for the vibrational energy: r = B1B2B3, and w = BppB. 



(6) 
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Sj^, determine the volume element in the collective space. The moments of inertia are calculated from the Inglis- 
Bclyacv formula: 
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where k denotes the axis of rotation, and the summation runs over the proton and neutron quasiparticle states. 
The mass parameters associated with the two quadrupole collective coordinates = {Q20) and 52 = {Q22) arc also 
calculated in the cranking approximation 
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Finally, the potential 14oii in the collective Hamiltonian Eq. (3) is obtained by subtracting the zero-point energy 
corrections from the total energy that corresponds to the solution of constrained RMF-I-BCS equations, at each point 
on the triaxial deformation plane. 

The Hamiltonian Eq. (3) describes quadrupole vibrations, rotations, and the coupling of these collective modes. The 
corresponding eigenvalue problem is solved using an expansion of cigenfunctions in terms of a complete set of basis 
functions that depend on the deformation variables /3 and 7, and the Eulcr angles (f>, 9 and [21]. The diagonalization 
of the Hamiltonian yields the excitation energies and collective wave functions: 



^i^\Pn,^)= E €Al3n)'^iiKm- (10) 



KeAi 

The angular part corresponds to linear combinations of Wigner functions 



and the summation in Eq. (10) is over the allowed set of the K values: 

A/=(2'''---'{ , (12) 

[2,4, ...,/- 1 for / mod 2 = 1 . ^ ' 

In Figs. 5 and 6 we plot the isotopic dependence of two characteristic collective obscrvablcs: i?4/2 and B(E2; 
2^ ^0]'") (in Weisskopf units), for Ba and Xe nuclei, respectively. The values obtained by the diagonalization of the 
collective Hamiltonian Eq. (3), with parameters determined by the microscopic relativistic energy density functional 
PC-Fl, are shown in comparison with data [24, 25]. For a transition between 7-soft rotors (0(6) dynamical symmetry 
limit of the interacting boson model) and spherical vibrators {U{5) dynamical symmetry limit), the ratio between the 
excitation energies of the first 4+ and 2+ states varies from the value R4/2 — 2.5 in the 0(6) limit, and R4/2 = 2 for a 
spherical vibrator. For a second-order shape phase transitions between these limits, the parameter-free E(5)-symmetry 
prediction at the critical point is R4/2 = 2.2. The calculation reproduces the rapid decrease of R4/2 with mass number. 
The agreement with the empirical ratios is somewhat better in the case of Ba isotopes, but even for Xe nuclei the 
differences are not large. The largest deviation from the experimental value is calculated for ^"^"Xe: « 0.15. For the 

= 78 isotones ^'^^Ba and ^^^Xe, the calculated ratios R4/2 ~ 2.3 are not very different from the E(5)-symmctry 
prediction. An even better agreement with data is obtained for the calculated isotopic dependence of the transition 



5 



probabilities B(E2; 2^ — !• 0'^) (panels on the right of Figs. 5 and 6). For instance, the calculation reproduces in 
detail the swift decrease of the empirical B(E2) values from « 40 Weisskopf units in ^^*Xe, to about 15 Weisskopf 
units in ^'^'^Xe. It should be noted that the calculation is completely parameter-free, i.e. physical observables, such 
as transition probabilities and spectroscopic quadrupole moments, are calculated in the full configuration space and 
there is no need for effective charges. 

Before comparing the calculated excitation spectra and E2 transition rates with available data and E(5)-symmetry 
model predictions, in Figs. 7 and 8 we plot the relative fluctuations of the quadrupole deformations A/3/ < (3 > and 
A7/ < 7 >, for the sequence of ground-states of Ba and Xe isotopes, respectively. The variances are defined by [26]: 



are calculated in the nuclear ground state, i.e. in the lowest eigenstate of the collective Hamiltonian. It is interesting 
to note the change in the fluctuations around N ~ 78. For a transition from deformed to spherical shapes one expects 
an increase of the relative fluctuation A/3/ < /3 >. For the triaxial deformation variable 7 the calculation predicts a 
marked maximum of the relative fluctuation in ^"^"^Ba and ^'^^Xc, that can be interpreted as a possible fingerprint of 
the second-order shape phase transition. 

^^^Ba was the first nucleus to be considered as a good example of empirical realization of the E(5) symmetry at 
the critical point of second-order phase transition between spherical and 7-soft shapes [8] . In Fig. 9 we compare the 
spectrum of the collective Hamiltonian for -'^"^'^Ba with available data for positive parity states [24, 25, 27, 28], and 
with the predictions of the E(5) model. For the moments of inertia of the collective Hamiltonian we have multiplied 
the Inglis-Belyaev values Eq. (7) with a common factor determined in such a way that the calculated energy of the 2f 
state coincides with the experimental value [19, 21]. This additional scale parameter is necessary because of the well 
known fact that the Inglis-Belyaev (IB) formula (7) predicts effective moments of inertia that are considerably smaller 
than empirical values. Following the minimal prescription of Ref. [29], the moments of inertia that parameterize the 
collective Hamiltonian can be simply related to the IB values: Xk{q) = ^k^{<l){^ + c^), where q denotes the generic 
deformation parameter, and a is a constant that can be determined in a comparison with data [19]. Below we discuss 
in more detail the approximations used in the calculation of moments of inertia and mass parameters. 

The excitation spectrum of the collective Hamiltonian determined by the PC-Fl density functional is in good 
agreement with the available data. This is especially true for the excitation energies and transitions in the ground- 
state band. The calculated sequence of states 2 J, 3^, and 4^, is approximately 300 keV higher than the corresponding 
experimental levels, and the transition 2^ — > 2'^ is considerably weaker compared to experiment (notice, however, 
the large error bar of the experimental value). By comparing the de-excitation pattern of the two excited 0"*" states, 
it appears that the ordering of these states is reversed in the theoretical spectrum. The calculation predicts that the 
second 0"*" state predominantly decays to 2f, and for the third 0+ state a much stronger transition is predicted to 2 J. 
In fact, the calculated Oj state is close in energy to the experimental 0|^ state which displays a similar E2 branching 
to 2i and 2 J. The third theoretical O"*" state, however, is calculated more than 1 MeV above the corresponding 
experimental state. One might also notice that the absolute B(E2) values (in Weisskopf units) for the stronger 
transitions are in very good agreement with data, but the calculated Oj and 0'^ states are obviously much more mixed 
than the corresponding experimental states. This also explains the strong repulsion between the theoretical second 
and third 0+ states. 

It is interesting to compare the PC-Fl excitation spectrum with the E(5) level scheme that corresponds to ^^''Ba 
(right panel in Fig. 9). In addition to the total angular momentum /, in the case of E(5) symmetry states are 
labeled by two quantum numbers: ^ = 1,2,3,... enumerates the zeros of the Bessel functions Jr+3/2 in /3, and 
T = 0,1,2,... is associated with the 0(5) algebra and denotes the phonon-like levels within a ^-family [7]. The 
corresponding spectrum is parameter-free up to an overall scale factor, that is adjusted by normalizing the energies 
to the experimental excitation energy of the first excited state 2f. The E(5) transition rates are normalized to the 
experimental B(E2; 2^ — )■ O]*"). Some striking similarities are found in the comparison of PC-Fl and E(5) excitation 
spectra for ^'^'*Ba. For instance, in the ground state band the calculated ratios 




(13) 



where the average values of /3 and 7 



< 13 >^ ^ (/32), < 7 >= arccos 
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= 1.67 
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B{E2;A+ 2+) 



1.33 , 



exactly coincide with the values predicted by the E(5) symmetry model. Again, we emphasize that the calculation 
of transition probabilities in the microscopic collective model is parameter-free. A very good agreement is also found 
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for B(E2; 4^ ^- 2^) and B(E2; 3^ ^ 2^). The second 0+ state is calculated at an excitation energy very close to the 
predicted position of 0^^2 r=0' ^^'^ state also predominantly decays to 2^ (At = 1 transition in the E(5) model). 
Even though calculated at much higher energy, the decay pattern of the third O"*" state is similar to that of the third 
E(5) 0+ state with ^ = 1, and r = 3. The two excited 0+ states in the E(5) model belong to different ^-multiplets, and 
only At = 1 transitions are allowed. The microscopic collective Hamiltonian, on the other hand, breaks the symmetry 
and this is reflected both in the decay pattern and the excitation energy of the O3 state. Another difference between 
the PC-Fl and E(5) spectrum is the transition 2^ -> 2f, for which the E(5) model predicts the same B(E2) value as 
for 4^ — >■ 2^. The smaller value calculated from the eigenvectors of the collective Hamiltonian can be understood by 
considering the corresponding probability density distributions. For a given collective state Eq. (10), the probability 
density distribution in the (/3, 7) plane is defined by 

P/a(/3,7)= E IV'aK(/3,7)P/3'|sin37| . (15) 
KeAi 

The normalization reads: 

/ /3d/3 / d7P/a(/3,7) = l- (16) 
Ja Jo 

In Fig. 10 we plot the probability densities in the /3 — 7 plane for the states 2f, 2^ and of -'^^^Ba. The states 2f 
and 4^ belong to the ground-state band, and the corresponding wave functions display a large overlap in the /3 — 7 
plane. The average values of (3 are: 0.156 for 2f, and 0.166 for 4'^, respectively, and the average values of 7 are: 
23.4° for 2^ , and 23.7° for 4^. Because of the large overlap of the wave functions, one expects a strong transition 
between these states. For the state 2^ the average value of of /? is 0.174, but < 7 >= 29.6°, i.e. the overlap between 
the collective wave functions of the states 2^ and 2^ is smaller than the one between 2^ and 4^. Consequently, the 
transition 2 J — > 2f is weaker than 4^ — >■ 2^. However, as we have already noted, this result is not in agreement with 
data which show that the transition 2 J — 2^ is in fact stronger than 4^ — )- 2]^. 

Turning now to the sequence of Xe isotopes, in Fig. 11 the PC-Fl collective excitation spectrum of ^^^Xe, i.e. the 
isotone of ^'^'*Ba, is compared with available data [24, 30], and the level scheme predicted by the E(5) symmetry 
model. We find reasonable agreement with experiment and, again, a remarkable similarity with the decay pattern 
predicted for the E(5) critical-point symmetry. The limited set of data, however, do not provide strong evidence for 
^^^Xe being a good example of empirical realization of the E(5) critical-point symmetry. Therefore in Fig. 12 we 
also display the spectra of ^^°Xe. In this case the microscopic calculation based on the PC-Fl functional predicts an 
excitation spectrum in slightly better agreement with available data [24, 30], e.g. the transitions between low- lying 
states, and the decay scheme is also very similar to the one that characterizes the E(5) symmetry, except for a weaker 
transition 2^ -> 2^ that has already been discussed in the case of ^'^^Ba. One might even notice that the calculated 
E2 transitions from Oj and 0'^ are, for this nucleus, actually closer to the branching pattern predicted by E(5) than 
they were in the case of ^'^"'Ba and ^^^Xe. The transitions 0;^ — >■ 2f and 0^ — >■ 2^, which are forbidden in the E(5) 
symmetry limit because of the Ar — 1 selection rule, are strongly suppressed between eigenstates of the PC-Fl-bascd 
collective Hamiltonian. 

In a systematic search for experimental fingerprints of E(5) critical-point symmetry [31], it was suggested that 
one possible candidate might be ^^®Xe. However, for this nucleus the ratio R4/2 = 2.33 is intermediate between the 
value for E(5) {R4/2 = 2.2) and the deformed 7-independent limit or 0(6) {R4/2 = 2.5). In fact, the microscopic 
potential of ^^^Xe (Fig. 2) displays a 7-independent surface for axial deformations 0.1 < f3 < 0.2. We also notice 
that the potential is not nearly as rigid in /3 for /3 < 0.1 and (3 > 0.2, as one assumes in the E(5) symmetry limit. 
In Ref. [32] a 7-soft analog of the confined /3-soft (CBS) rotor model [33] has been constructed, that corresponds to 
an infinite square-well potential over a confined range of values > Z?™ > (confined /3-soft potential, jSm and 
Pm denote the positions of the inner and outer walls on the (3 axis, respectively). The 7-soft extension of Ref. [32] 
generalizes the E(5) critical-point symmetry to a parametric solution over the whole path between E(5) and the /3-rigid 
deformed 7-independent limit. The structure of ^^*Xe was investigated and, from the decay pattern of the first two 
excited 0+ states, it was suggested that ^^*Xe is located well beyond the U{5) — ^ 0(6) shape-phase transition in the 
deformed phase. In a very recent experimental study of low-lying collective states in ^^*Xe by 7-ray spectroscopy 
[34], the data, including the first measurement of the B(E2) values for transitions from the Oj and 0;^ states, were 
compared with theoretical predictions and it was concluded that ^^^Xe is not an E(5) nucleus, thus leaving ^•^°Xe as 
the best candidates among the Xe isotopes. In Fig. 13 we, therefore, finally compare our microscopic results for ^^*Xe 
with the recent experimental decay scheme of Ref. [34]. The calculated excitation energies and B(E2) values for the 
ground-state band and the sequence 2^, 3^, 4^, are in good agreement with experiment. The two excited 0"*" states, 
however, are calculated at higher energy, especially 0;^. The predicted transitions to 2f and 2^ clearly follow the 
pattern calculated in -'^'^^Xe and ^^°Xe, and do not reproduce the empirical two-level mixing that strongly breaks the 
At = 1 selection rule characteristic for 0(5) symmetry. 
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For all nuclei analyzed in this work we have found a very good agreement between data and the predictions of 
the microscopic collective Hamiltonian based on the PC-Fl energy density functional plus (5-force pairing, except for 
the two excited 0+ states. In particular, the second excited 0+ state, i.e. the one that predominantly decays to 2^, 
is systematically calculated at considerably higher excitation energy compared to the experimental position of the 
two excited 0"*" levels. In ^^"'Ba and ^^®Xe, for which data on absolute E2 transition strengths are available, it is the 
first excited 0+ state that predominantly decays to 2 J, contrary to the results of our calculation. One reason for 
this discrepancy, of course, could be the particular choice of the energy density functional and/or the treatment of 
pairing correlations, but also the approximation used in the calculation of the mass parameters (vibrational inertial 
functions). In the current version of the model the moments of inertia are calculated from the Inglis-Belyaev formula 
(Eq. (7)). These moments of inertia do not include the Thouless-Valatin dynamical rearrangement contributions 
and, therefore, are systematically smaller than the empirical values. However, as shown in several studies [29], the 
Thouless-Valatin corrections are almost independent of deformation, and the effective moments of inertia to be used 
in the collective Hamiltonian can simply be obtained by renormalizing the Inglis-Belyaev values by a constant factor, 
characteristic for a given nucleus. The situation is considerably more complicated in the case of mass parameters 
[35, 36], for which there are no simple estimates of the Thouless-Valatin correction, especially for nuclei with 7-soft 
potential energy surfaces. Some authors [36] argue that, to approximately take into account the Thouless-Valatin 
correction, all inertial functions, not only the moments of inertia, should be rcscaled by a constant multiplicative 
factor. 

Another effect that is not included in our model is the coupling of nuclear shape oscillations with pairing vibrations, 
i.e. vibrations of the pairing density. A number of studies have shown that excited 0+ states are very sensitive to 
the coupling between these modes. For the Ba and Xe isotopes, in particular, the effect of coupling between pairing 
vibrations and axial quadrupole vibrations was investigated in Ref. [37], using collective Hamiltonians based on a 
microscopic Nilsson single-particle Hamiltonian with pairing forces and a long-range residual interaction. It was 
shown that the coupling between the pairing and quadrupole modes has a very pronounced effect on the lowest 
excited 0^ states, lowering their excitation energies by 40 — 50%, in much closer agreement to data. In an even earlier 
study of collective states in Xe and Ba isotopes [38], based on the general 5-dimensional Bohr Hamiltonian, this 
coupling was not taken into account explicitly, but rather it was simulated by a reduction of the pairing strengths of 

20%. Even though the resulting pairing interaction is somewhat unphysical, i.e. it gives proton and neutron gaps 
that are too small, the resulting excitation energies of the lowest 0^ states are 40 — 50% lower, in better agreement 
with experiment. 

The present study, however, employs exactly the same model that was used in our previous analysis of spherical 
to axially-dcformed shape transitions, The mass parameters are calculated in the cranking approximation Eqs. (8) 
and (9), and from these expressions it is obvious that the resulting mass parameters will be very sensitive to the 
underlying microscopic shell structure and to the treatment of pairing correlations. In fact, we have found that the 
position of the first excited 0+ state in our calculation (the one that predominantly decays to 2^) is sensitive to the 
value of the mass parameter Bpfj, the excitation energy of the third 0^ state (predominantly decaying to 2 J) displays 
a very strong dependence on B^^, and the mixing of the two excited 0"*" states is controlled by Bp^. For instance, by 
decreasing B/j^^ i.e. by reducing the coupling between /3 and 7 vibrational degrees of freedom, the mixing between 
O2 and 0;^ decreases and the resulting spectrum is even closer to the E(5) level scheme. By increasing the value of 
Bj^ one finds a rapid decrease of the excited energy of the calculated O3 , and eventually the two excited 0+ state 
cross, in better agreement with available data for ^'^^Ba and ^^^Xe. The only other effect on the calculated spectrum 
is a much smaller lowering of the sequence: 2^, 3^, 4^, again in agreement with experiment. In the case of ^^^Xe, for 
instance, an enhancement of B-y-y by a factor > 1.5 brings the calculated excitation energies in very good agreement 
with the experimental spectrum shown in Fig. 13 [34]. 

Very recently a global study of low-energy nuclear structure at normal deformation has been reported, based on 
the non-relativistic Hartree-Fock-Bogoliubov framework with the Gogny interaction DIS, and mapped onto a five- 
dimensional collective quadrupole Hamiltonian [39]. Ground-state properties, excitation energies, quadrupole and 
monopole transition matrix elements for yrast levels up to 6^ and the lowest yrare levels, have been calculated for 
even-even particle-stable nuclei with proton numbers Z = 10toZ=:110 and neutron number TV < 200. Except for 
the use of a different effective interaction in the particle-hole and particle-particle channels, the model of Ref. [39] is 
very similar to the one employed in the present study. However, while in both models the quadrupole mass parameters 
are calculated in the cranking approximations, in the calculation of the moment of inertia the model of Ref. [39] goes 
beyond the simple Inglis-Belyaev formula used in this study, and explicitly includes the Thouless-Valatin dynamical 
rearrangement contributions. This leads to an increase of the moments of inertia and the resulting excitation energies 
of yrast states in deformed and transitional nuclei are, in general, in very good agreement with data. As a particular 
example, the study of Ref. [39] has considered the low-energy collective quadrupole spectrum of ^^^Sm. It has been 
shown that the spectrum calculated with the collective Hamiltonian based on the Gogny interaction DIS is consistent 
with the X(5) analytical solution at the critical point of QPT between spherical and axially deformed shapes in 
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iV 90 isotones. More generally, the calculated spectral properties of « 1700 even-even nuclei have demonstrated 
the predictive power of the quadrupole collective Hamiltonian model based on a microscopic effective interaction. A 
problem common to both studies, i.e. to the global survey of Rcf. [39] and the specific calculation of transitional Ba 
and Xe isotopes reported in this work, are the excitation energies of excited 0+ states that are systematically calculated 
too high compared with data. Similar to the discussion above, in Ref. [39] it has been suggested that the energies 
of excited 0"^ states could be improved by including the Thouless-Valatin dynamical rearrangement contributions in 
the calculation of the quadrupole mass parameters, and by extending the model space to include two-quasiparticle 
components in the wave functions. 

III. CONCLUSIONS 

The recently developed implementation for the solution of the five-dimensional collective Hamiltonian for quadrupole 
vibrational and rotational degrees of freedom, with parameters determined by constrained self-consistent relativistic 
mean-field calculations for triaxial shapes, has been applied to the study of 7-soft shapes in Ba and Xe nuclei in the 
mass A > 130 region. This microscopic approach had previously been used to describe spherical to axially deformed 
shape transitions in the mass A w 150 region, and the results have shown there is an abrupt change of structure at 

= 90 that can approximately be characterized by the X(5) analytic solution at the critical point of a first-order 
quantum phase transition. In this sense the present analysis represents a natural extension of our previous studies 
[15, 19, 20], to a region of nuclei where evidence has been reported for a second-order QPT between spherical and 
7-soft shapes, with the critical point of phase transition characterized by the E(5) dynamical symmetry. Without any 
modification of the model Hamiltonian, of the treatment of pairing correlations, or the parameters of the effective 
interactions in the particle-hole and particle-particle channels, we have been able to describe the rapid transition 
between spherical and 7-soft shapes in Ba and Xe isotopes. The results reproduce the characteristic evolution of 
excitation spectra and E2 transition probabilities, and a good agreement with available data is obtained. In specific 
nuclei, e.g. ^■^^Ba and ^■^^Xe, the calculated triaxial quadrupole binding energy surfaces are rather flat in an extended 
interval of the axial deformation parameter /3, and almost independent of the 7-deformation. For the ground states of 
these nuclei one finds pronounced discontinuities in the relative fiuctuations of the quadrupole deformations A/3/ < /3 > 
and A7/ < 7 >. The calculated spectra are remarkably similar to the level schemes predicted by the E(5) dynamical 
symmetry model and, therefore, display fingerprints of a second-order shape phase transition. 

The fact that the same microscopic model describes both first- and second-order QPT in different mass regions 
is, in our opinion, an interesting result which illustrates the importance of microscopic studies of QPT starting from 
nucleonic degrees of freedom. A quantitative comparison with experimental spectra, particularly with very recent 
data for ^^®Xe, points to some intrinsic problems in the description of excited 0+ states and clearly indicates that one 
must go beyond the simple cranking approximation for the model mass parameters (vibrational inertial functions), 
or even extend the model by explicitly taking into account the coupling of nuclear shape oscillations with pairing 
vibrations. 
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FIG. 1. (Color online) Self-consistent RMF+BCS triaxial quadrupole binding energy maps of the even-even ^'^®~^^°Ba isotopes 
in the /3 — 7 plane (0 < 7 < 60"). All energies are normalized with respect to the binding energy of the absolute minimum, 
the contours join points on the surface with the same energy (in MeV). For each nucleus in the corresponding inset we plot the 
axial projection of the binding energy for oblate (negative (3) and prolate (positive /3) deformations. 
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FIG. 3. (Color online) Self-consistent RMF+BCS binding energy curves of the ^^''Ba nucleus, as functions of the deformation 
parameter 7, for four values of the axial deformation /3 = 0.05, 0.1, 0.15, and 0.2. 
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FIG. 4. (Color online) Same as described in the caption to Fig. 3 but for the nucleus ^'^^Xe. 




FIG. 5. (Color online) Evolution of the characteristic collective observables i?4/2 and B(E2; 2f Of) (in Weisskopf units) 
with mass number in Ba isotopes. The microscopic values calculated with the PC-Fl energy density functional are shown in 
comparison with data [24, 25]. 




FIG. 7. (Color online) Fluctuations of the quadrupole deformations A/3/ < P > and A7/ < 7 > for the sequence of ground- 
states of Ba isotopes. 
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FIG. 8. (Color online) Same as described in the caption to Fig. 7 but for the isotopes of Xe. 
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FIG. 9. (Color online) The low-energy spectrum of ^''^Ba calculated with the PC-Fl relativistic density functional (middle), 
compared with the data (left), and the E(5)-symmetry predictions (right) for the excitation energies, intraband and interband 
B(E2) values (in Weisskopf units). The theoretical spectra are normalized to the experimental energy of the state 2f , and the 
E(5) transition strengths are normalized to the experimental B(E2; 2f — > Oj'"). 



19 




FIG. 10. (Color online) Probability densities Eq. (15) in the /? — 7 plane for the states 2]*", 2j and 4^ of ^''''Ba. 
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FIG. 13. (Color online) The experimental low-energy spectrum of ^^*Xe (left), compared with the level scheme and decay 
pattern predicted by the solution of the microscopic collective Hamiltonian with the PC-Fl relativistic density functional. 



